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Abstract 

A rational number can be naturally presented by an arithmetic computation 
(AC): a sequence of elementary arithmetic operations starting from a fixed con¬ 
stant, say 1. The asymptotic complexity issues of such a representation are studied 
e.g. in [2, 9] in the framework of the algebraic complexity theory over arbitrary 
field. 

Here we study a related problem of the complexity of performing arithmetic 
operations and computing elementary predicates, e.g. “=” or on rational 
numbers given by AC. 

In the first place, we prove that AC can be efficiently simulated by the exact 
semidefinite programming (SDP). 

Secondly, we give a BPP-algorithm for the equality predicate. 

Thirdly, we put ^-predicate into the complexity class PSPACE. 

We conjecture that ^-predicate is hard to compute. This conjecture, if true, 
would clarify the complexity status of the exact SDP — a well known open prob¬ 
lem in the field of mathematical programming. 

Keywords: semidefinite programming, complexity, succinct representation. 

Algorithmic complexity provides a general framework to analyze complexity of 
computational problems. It works for many cases and gives results that are impor¬ 
tant for practical applications. Nevertheless, some basic assumptions of the theory are 
strange from practical point of view. We are sure that a linear time algorithm running 
in time 10 1QOO n cannot be realized in our Universe. Also, some widely used algorithms 
have exponential running time in the worst case. 

Numerical algorithms are especially in striking disagreement with complexity the¬ 
ory. Probably, the most popular exponential algorithm is the simplex algorithm for 
linear programming. It is widely used despite the existence of polynomial algorithms 
that were found after pioneering breakthrough of L.Khachiyan (see, e.g. [8, 7, 13]). 

The complexity analysis of the semidefinite programming problem involves even 
more difficulties. SDP is often considered as tractable due to various approximate 
algorithms. SDP is a convex optimization problem so the ellipsoid method can be 
applied to solve it approximately as well as a variety of interior points methods [6, 10, 
14, 18]. But there are amazingly few results on complexity of the SDP problem. 
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Khachiyan and Porkolab [11] found polynomial time algorithm for the SDP prob¬ 
lem when the dimension is fixed. They established doubly exponential bounds on the 
solutions for the general SDP problem and on the discrepancies of infeasible programs. 

For our further considerations the most interesting is Ramana’s result [12], He de¬ 
veloped the exact duality theory for SDP. Ramana’s dual program can be constructed 
in polynomial time. His analogue of Farkas lemma has an immediate complexity- 
theoretic corollary: a complement to the SDP feasibility problem (SDFP) can be re¬ 
duced to SDFP itself. It implies that SDFP cannot be NP-complete unless NP = 
coNP. More generally, SDFP should belong to a complexity class that is closed un¬ 
der complement. Examples are P, NP fl coNP, BPP, PSPACE. (For definitions of 
these classes and other useful information on complexity theory see, e.g., the Sipser’s 
book [16].) 

Here we address the exact SDFP problem. The problem is to check that intersection 
of the cone of positive semidefinite matrices with some affine subspace of matrices 
is not empty. The subspace involved is defined by generators that are matrices with 
rational entries. It is well known that some pathological examples exist for the SDFP. 
For instance, it is possible that feasible program has only doubly exponential solutions 
and that an infeasible program could have doubly exponentially small discrepancy. 
These examples show that in some cases a polynomially bounded machine can not 
simply write a solution to SDFP. Of course, this does not prevent checking the existence 
of a solution in polynomial time. 

Motivated by these examples, we introduce a new problem: comparison of arith¬ 
metic computations. It seems to be interesting in it’s own right and concerns a nonstan¬ 
dard way of representing integers and rationals. The common way to represent integers 
is to use a positional system. Among positional systems the binary system is the sim¬ 
plest and the most natural. Binary representation of a number TV is a string of length 
6»(log TV). This bound is optimal due to a counting argument. But it is also possible to 
encode numbers in such a way that some numbers are encoded by very short strings. 
In this case we speak of a succinct representation of an integer. A natural way for suc¬ 
cinct representation of a rational number r is to use an arithmetic computation or an 
arithmetic circuit (AC). By definition, AC is a sequence of the elementary arithmetic 
operations starting from a fixed constant, say 1, and generating r as output. 

To our knowledge, the complexity issues of the AC representation were studied 
primarily in the asymptotic setting in the framework of the algebraic complexity theory 
over a field (see, e.g. [15, 9]). Namely, let r(k) be the minimum number of arithmetic 
operations required to build the integer k from the constant, say 1, i.e. the length of 
the minimal AC computing k. A sequence Xk of integers is said to be “ultimately 
easy to compute” if there exists another sequence ak and polynomial p(-) such that 
T~(akXk) ^ p(logfc) for all k (for instance, in can be shown that the sequence 2 k is 
ultimately easy to compute). Otherwise the sequence is said to be “ultimately hard 
to compute”. Counting argument shows that ultimately hard to compute sequences do 
exist. A central open problem is to show that some explicit sequences, say, n\, |_(3/2)”J 
are ultimately hard to compute. For instance, if n! is ultimately hard to compute then 
Pc = NPc over the field of complex numbers [15]. 

It is worth to mention that the complexity of AC is related to a circuit model of 
computation for polynomials that was introduced and studied by Valiant [17]. The 
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Valiant’s model turns to AC by substitution constants in polynomials. Some interesting 
applications of this observation are discussed in [9]. 

In this paper we explicitly treat a related problem of the complexity of performing 
arithmetic operations and computing elementary predicates, e.g. “=” or “Jj”, on ratio¬ 
nal numbers represented by AC. An immediate inspection shows that if numbers are 
repesented by AC then to perform an elementary arithmetic operation one should sim¬ 
ply merge the AC of the operands in the appropriate way, thus arithmetical operations 
are easy to perform, despite the evident fact that AC representation is succinct for some 
numbers (e.g. 2 2 can be obtained by repeated squaring). On the other hand, — it is 
unclear how to efficiently compute elementary predicates “=” or i.e. how to check 
efficiently whether two AC compute the same value or how to compute the maximal 
value. In the sequel we denote these predicates by AC = and ACj» respectively. 

Our main contributions are: we prove that some restricted version of AC can be 
efficiently simulated by the exact semidefinite programming (SDP) and construct a 
polynomial reduction of AC^ to SDFP; we give a BPP-algorithm for the equality 
predicate; we put ^-predicate into the complexity class PSPACE. 

Unfortunately, we are unable at the moment to prove any lower bounds for the 
^-predicate. Any nontrivial result in this direction would imply lower bounds for the 
exact SDP — one of the main open problems in the field of mathematical programming. 

The rest of paper is organized as follows. In Section 1 we give a definition of the 
circuit representation for rationals and state the basic results concerning it. Also, this 
Section contains technical results about a choice of a basis for AC. In Section 2 we 
construct a reduction of AC^ to SDFP. Sections 3, 4 contain proofs of AC = £ BPP 
and AC^ £ PSPACE. In the last Section 5 we discuss open questions around the 
AC^ and SDFP problems. 

1 AC representation of numbers 

Let B be a finite collection of functions of type Q fc —> Q (k may vary). It is called a 
basis. A circuit over basis B is a sequence S = so, Si,..., se of assignments such that 
so := 1 and for each i ^ 1 there exist j,k < i such that s, := Sj * Sk and * £ B. The 
size £(S) of circuit S is the number t. 

The most natural basis consists of four arithmetic operations {+, —, /, •}. Circuits 
over this basis are called arithmetic circuits. Circuits over the basis {+, —, •} are called 
division-free circuits. Monotone (arithmetic) circuits are circuits over the basis {+, •}. 

We will represent rationals by circuits. For each circuit we define the value of the 
circuit by induction. We will use notation v(S) for the value of circuit S. Value of 1 
is 1. The value of a circuit S = so, si, ..., se where se = Sj * Sk is v(Sj) * v(Sk). 
Here Sj = so, si,..., Sj. Note that each prefix of a circuit is a circuit by definition. If 
some operations can not be performed (e.g. a division by 0) the value of such a circuit 
is undefined. Let X = v(S). We say that such AC represents X. 

It is easy to see that arithmetic circuit representation can be much more compressed 
than the usual binary representation. A circuit 5 = 1, si, S 2 , ■ ■ ■ ,se where si = 1 + 1 
and Sj + 1 = Sj ■ Sj for j f 1 represents 2 2 . On the other hand this example is 

asymptotically optimal. The following statement can be easily verified by induction. 
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Statement 1. If p/q = v(S) where p, q are integers then maxjp, q} = 0( 2 2 ). 

1.1 The complexity of equality and inequality predicates over AC 

Implementation of arithmetic operations with rationals represented by arithmetic cir¬ 
cuits is straightforward and can be done in linear time. Formally we write 

<S(X * Y) = S(X), tail(<S(y)), S£(s(.y))+*(S(y)-i)> 
where s^s{x))+e(S(Y))-i := s^s(x)) * $t(S(x))+t(S(Y))~ 2- 

Here tail(«S) means a circuit S without the starting 1. 

But how difficult may be the computation of the equality predicate and of the in¬ 
equality predicate? We state these algorithmic problems formally. We always assume 
that a circuit is represented by a list of triples (*, j, k ) where * £ B and j, k are positive 
integers. The ?rth element in the list corresponds to an assignment s n = Sj * Sfc. Thus 
a circuit of size £ is written as a 0{l logf) binary word. 

Predicate AC = (B). It is true for a pair of circuits <S], S 2 over the basis B iff 

u(Si) = v{Sf). 

Predicate AC^ ( 13 ). It is true for a pair of circuits Si, S 2 over the basis B iff 
u(5i) > v(S 2 ). 

Note. If B contains division then these predicates are partially defined (value of 
some circuits can be undefined). Partially defined predicates are called promise prob¬ 
lems. Most complexity classes can be easily redefined to include promise problems 
and most results remain true in this, more general, setting. We omit a discussion of 
promise problems but indicate that their use is safe in our considerations. 

We denote the equality and the inequality predicates over the arithmetic basis by 
AC= (resp. AC^). 

It is clear that both predicates fall into the complexity class EXPTIME. In fact, 
we are able to place them into lower levels of computational hierarchy. 

Theorem 1. ACL G BPP. 

In other words, the equality check can be performed by probabilistic Turing ma¬ 
chine in polynomial time. 

We are unable to give an exact characterization of the computational complexity of 
the second predicate. Computing AC^ looks as a computationally hard problem. An 
obvious way to solve it is to make all calculations indicated in the circuits that form 
input of the problem. Using binary representation we need exponentially large memory 
to do it. Using modular arithmetic it is possible to solve AC^ in polynomial memory. 

Theorem 2. AC^ G PSPACE. 

The proof of Theorem 2 uses some constructions from an NC 1 algorithm for com¬ 
parison integers in modular arithmetic [4, 5], 
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1.2 Equivalent bases 

Two bases B\ B 2 are called “= ” equivalent (respectively, “M’ equivalent) iff the pred¬ 
icates AC = (£>i) and AC = (£> 2 ) (respectively, AC^(£>i) and AC^(£> 2 )) are mutually 
polynomially reducible. 

Theorem 3. The following bases are “=” and “f ” equivalent: arithmetic, division- 
free, monotone and {+, x i—> x 2 /2\. 

The last basis in the list is added for technical purposes. It is used in the reduction 
of the problem ACj> to the feasibility problem for semidefinite programming. 

Reductions of all mentioned bases to the arithmetic basis are straightforward. To 
prove Theorem 3 we establish reductions in the opposite direction. 

Note that if a basis contains a subtraction then general predicate AC^(B) is re¬ 
ducible to its particular case when one of the compared values is zero. Indeed, suppose 
we are going to compare and v(S 2 ). We can merge the circuits Si, S 2 into one 

circuit. This merged circuit contains all assignments of Si, S 2 and ends by the assign¬ 
ment d := a — b, where a and b are the last assignments in the circuits <Si, £ 2 - The 
same argument can be also applied to the predicate AC—(B). 

All reductions described below have similar form. A circuit S over some basis is 
converted to a circuit S' over another basis using step-by-step substitution of constant¬ 
sized groups of assignments instead of each assignment in S. 

Lemma 1. AC^({+, —, /, •}) (resp. AC = ({+, —, /, •})) is reducible to AC^({+, —, •}) 
(resp. AC^({+, •})). 

Proof. Informally speaking, the lemma is very simple: we can keep and transform 
numerators and denominators separately. Below we present a more detailed description 
of the reduction. 

Let S be a circuit of size t over the arithmetic basis. We construct a circuit S' of 
size 0{t) over division-free basis in the following way. 

The circuit S' consists of four sequences of assignments A 1 , A 2 , M, V. For the 
assignment .s, := Sj ± Sk in S add the assignments 

A\ '■= A fj -T>k, A 2 := Vj • A4, T^i := T^j -Vk, M := A] ±A 2 . 

Similarly, for the assignment .s, := s 3 ■ Sk in S add the assignments 

T>i := Vj • V kl Mi := Mj ■ M k , 

and for the assignment := Sj/s k in S add the assignments 

Vi := Vj ■ M k , Mi := Mj ■ V k . 

The last assignment in the circuit S' is 

s'n ■= Me ■ Ve¬ 


it is easy to see that 


v{S) 


v(Me) 
v(V ( )' 
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So,v(S) ^ Oiffs^ > 0. 

Note that due to our assumptions the case Vi = 0 is impossible. So, the same 
reduction is valid for the predicate AC = . □ 

Lemma 2. AC^({+,— ,•}) (resp. AC = ({+, —, •})) is reducible to AC^({+,-}) 
(resp. AC^({+, •})). 

Proof. Informally, we do the same trick as above using equalities 

(. A-B) + (C-D) = (A + C)-(B + D ), (1) 

(. A-B)-(C-D) = (A + D)-(B + C ), (2) 

(A — B) ■ (C — D) = (A • C + B • D) — (B • C + A ■ D). (3) 

Now let consider the details of the reduction. 

At first we convert an input (<S l5 S 2 ) of ACj>({+, —, /, •}) into (<S, 0) as it was 
explained above. 

Then we construct two circuits £, R such that v(S) = v(C) — v(R). So, v(S) f 0 
iff v(C) ^ v(R) as well as v(S) = 0 iff v(C) = v(R). Again, the same reduction 
will work for both predicates. The size of £, R will be 0(s(S)) and they will be the 
circuits over {+,•}. 

Both circuits C and 'R. consist of six sequences of assignments C 1 , C 1 , C\ 7 Z 1 , R 2 , 
R?. 

For the assignment s-i := Sj + Sk in S add assignments to 1Z 

l\ := c) + 4 , n\ -.= n)+ n\ 

(see Eq. (1)). 

For the assignment s, := -Sj — Sk in S add the assignments to R 

4 :=C]+Rl R\ :=4+4 

(see Eq. (2)). 

For the assignment S{ := -Sj ■ Sk in S add the assignments to R 

4 := 4 • 4 > 4 : = n ) ■ 4 : = 4 + 4 . 

R-i ■■= n) ■ 4 . Ri ■■= 4 • Kl R\ := R- + R\. 

(see Eq. (3)). 

The circuit £ has the same structure except £, and R, assignments (or groups of 
assignments) are interchanged. By induction, we see that for each i the value of .s, is 

44) -v(Rj). □ 

Lemma 3. AC^({+, •}) (resp. AC=({+, •})) is reducible to AC^({+,a: x 2 /2}) 

(resp. AC = ({+,a: 1 —> a; 2 /2})). 

Proof. Let S be a circuit over the monotone basis. We construct a circuit S' over the 
basis {+,x 1 —> ar/2} consisting of 10 series of assignments V, Af, A 1 , t € [1,8] such 
that for each k the equality 

v(sk) = v(Vk) - v(Nk) (4) 
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holds. The construction is step-by-step substitution as in the above lemmas. 

For the assignment s, := Sj + s k in S add the assignments 

Vi := Vj + Vk i Af-i := Afj + Af k 

and for the assignment S{ := Sj ■ s k add the assignments 

A\ := Vj + Vk, A 2 := A/} + A 4, A® := Pj + A4, A^ : A/} + Vk, 

Al := {A\f/ 2, 4 := (A 2 ) 2 /2, Aj := (A®) 2 /2, Af := (^) 2 /2, 

Vi = A 5 i+ Al 
A U = Aj + A®. 

In the latter case the following equations hold 

v(Vi) = (v(Pj) + v(V k )) 2 /2 + (v(Afj ) + f(A4)) 2 /2, 
v(Afi) = (v(Vj) + v{Af k )f/2 + (w(A/}) + v(V k )f/2. 

Eq. (4) is verified by induction using Eq. (7) 

v(si ) = (v{Vj) - v(Afj))(v(V k ) - v(Af k )) = 

(v(Vj) ■ v{Vk) + v(Afj) ■ v(Af k ))— 

(v{Vj) ■ v(Af k ) + v(Afj) ■ v(V k )) = 

{(v{V j )+v{V k )f/2 + {v{Afj)+v{U k )f/2)- 
(( v (Vj) + v{Af k )) 2 /2 + (v{Afj) + v(V k )) 2 /2) = 

v(Vi)-v{Afi). (8) 

Now we are able to construct a reduction. Take an instance (Si, S?) of the prob¬ 
lem AC^({+, •}). Convert the circuits <Si, S -2 into the circuits iS 1 , S 2 over the basis 
{+. x i— > ;e 2 /2 } as described above. Join them into one circuit S. The circuit S' is 
an extension of S by the assignment f := V l Sl + Af 2 2 . Similarly, the circuit S" is an 
extension of S by the assignment f" := V 2 2 + A/"* . Here si is a size of <Si and S 2 is a 
size of S- 2 - 

A reduction is given by the mapping (<S-i. ^ 2 ) 1 —> (S', S”). From Eq. (4) we see 
that it also works for both types of predicates mentioned in the lemma. □ 

2 Reduction of the problem AC^ to SDFP 

Linear optimization on the intersection of the cone of positive semidefinite matrices 
with an affine subspace of matrices is called semidefinite programming (SDP). Let de¬ 
note SDFP the corresponding feasibility problem: to check whether the cone of positive 
semidefinite matrices has nonempty intersection with affine subspace of subspace. 
Semidefinite feasibility problem (SDFP) can be stated in the following form. 
Input: a list Qq, ..., () m of symmetric (n x n) matrices with rational entries. Ma¬ 
trices are represented by lists of entries, each entry is represented by a pair (numerator, 
denominator), integers are given in binary. 


(5) 


(6) 


(7) 
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Output: ‘yes’ if there exist reals x \,..., x m such that Q = Qo + YL'iL i x iQi is a 
positive semidefinite matrix, otherwise the output is ‘no’. 

The proof below uses Ramana’s results on the exact duality theory for SDP [12], 
Ramana found a special form of dual program for SDR It is called the extended La¬ 
grange - Slater dual program (ELSD). Ramana proved that 

• ELSD can be constructed from the primal program in polynomial time; 

• if both the primal and the dual are feasible then their optimum values are equal. 

Theorem 4 . The problem AC^ ({+, x h-> x 2 /2 }) is reducible to SDFP. 

Applying Theorem 3 we have reductions of the problem ACj> over all bases dis¬ 
cussed in the previous section to SDFP. 

Proof. The first step is to represent a circuit value as an optimal value of some semidef¬ 
inite program in the form 

t —> inf, 

(9) 

semidefinte conditions on t and other variables. 


Let S be a circuit over the basis {+,£ i—> x 2 /2}. We construct a semidefinite 
program P{S) such that for each assignment s; there is a variable Xi in P(S) and a 
matrix M,;. Matrices Mj are built as follows 

Si : = Sj -f- Sk > Mi — (Xi Xj X&) 

sr.= ,V2 - m ,= ( 2 « *>) 

Note that in the former case Mj >; 0 iff Xi ^ Xj + Xk and in the latter Mi A 0 iff 

Xi ^ x 2 /2. 

Let £ be a size of S. Consider a SDP 

X( —>• inf, 

Mi A 0 for each i, (10) 

x 0 = 1. 


To convert the program (10) to the standard form we replace an equality xq 
inequality 


(xo- 1 

V o 


0 ) 

— Xq + 1 J 


h 0 . 


1 by 

( 11 ) 


The optimal value of P(S) is the value of S. Indeed, all operations in the process of 
circuit evaluation are monotone with respect to each variable. So, any feasible solution 
of (10) satisfies the condition xe ^ v(S). On the other hand, Xi = v(si) is a feasible 
solution. 

Take now an instance of AC^({+,x i—> a: 2 /2}). It’s input is a pair of circuits 
5 (2 ). The positive answer in the problem AC^({+,x i—> at 2 /2}) means that 



^ u(<Sl 2 )) which is equivalent to —v(S^) ^ —u(5^ 1 - ) ). The value — v(S^) 
is an optimal value for SDP 

( 2 ) 

— x\ -» sup, 

M} 2) A 0, (12) 

4 2) = !• 

To represent -v(S rl 1 ) as the infimum of SDP we use the extended Lagrange- 
Slater dual program to the program of type (10). In Ramana’s paper ELSD has a mixed 
form (linear equations are permitted). To convert it to the standard form each linear 
equation should be replaced by a positive semidefinite condition on a (2 x 2) matrix as 
in Eq. (11). After convertion ELSD can be written as follows 


c(Y) —»inf, 
DAO. 

Here c(Y) is a linear functional on dual variables Y. 
Thus, the program 


-xf^ - c(Y ) > 0, 

Aff } A 0, 

x (2) - 1 
DAO. 

is feasible iff ^ 

This gives a reduction of AC^({+, x i—> x 1 j2\) to SDFP. 


(13) 


(14) 


□ 


3 The proof of Theorem 1 

The idea is to check the equality modulo random number. 

It is easy to compute a remainder of a circuit value modulo ‘short’ number (repre¬ 
sented in binary). Addition, multiplication and integer division are made in polynomial 
time by standard algorithms. 

We show that if two circuits have different values then with big enough probability 
they have different residues modulo random number. 

We present a BPP-algorithm for the monotone basis. By theorem 3 we conclude 
that AC = € BPP. 

The probabilistic algorithm for AC = ({+, •}). 

Input: a pair («Si, S 2 ) where Si, S> are circuits over the monotone basis. Let £ be the 
maximal size of the circuits. 

Step 1. Set B to 2 2 ^. 
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Step 2. Choose a random integer m from the uniform distribution on the interval 

Step 3. Compute r± = u(5i) mod m, r 2 = u(«S 2 ) mod m by making all operations 
indicated in the circuits <S|, iS 2 modulo M. 

Step 4. If r\ = r 2 then output ‘u(<Si) = n(<S 2 )’ else output ‘v(S\) ^ v[Sf)' ■ 

It is clear from the description that the algorithm uses 2£ random bits and runs in 
polynomial time. 

Claim 1. If u(iSi) = v(<S 2 ) then the algorithm outputs ‘u(iSi) ^ n(<S 2 )’ with 
probability 0. 

This is clear. 

Claim 2. If u(<Si) ^ u(<S 2 ) then the algorithm outputs ‘u(iSi) ^ i>( 5 2 )’ with 
probability at least (2£) _1 . 

Claim 2 is derived from Lemma 4 stated below. 

To fit the common definition of the class BPP we need to amplify success proba¬ 
bility by the standard procedure of multiple repetitions of the algorithm. Claims show 
that the gap to be amplified is I2(n _1 ) where n is an input size. So, an amplification can 
be made by poly(n) repetitions and the resulting algorithm runs in polynomial time. 

Now we need a lower bound of the least common multiple of integers taken from 
an interval [1, B] provided that we take sufficiently many integers. 

Lemma 4. Let 

N e {B) = min{lcm(xi,..., x r ) : Xi € [1 ,B],r> (1 — e)B}. 

Then for all s < (21nf3) _1 we have 

N e (B) > 2 n ( B / lnB ). 

Proof. Consider a set of integers R = {x \,..., x r } C [1, B\ such that |f?| > (1 —e)B. 
Factorize the least common multiple of these integers 

X = lcm(a;i,..., x r ) = p^ 1 ■ ■ ■ p a s ‘ ■ 

We are going to show that s = 12 (B/ In B). Then the lemma will follow from trivial 
bound X > 2 s . 

The prime number theorem gives an asymptotic 

Tl 

tt( n ) ~ - (15) 

Inn 

for the prime counting function u(n) = #{p : p ^ n,p prime}. 

Let R = [1, B] \ {xi ,..., x r }. At least ir{B) — s primes in the interval [1, B\ 
belong to the set R. From an inequality tt( B) s f B r we conclude that 

(1 — e)B < r < B — tt(B) + s 

and s = Q(B/lnB) fore < (2\nB)~ 1 . □ 
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To complete the proof of Theorem 1 we derive Claim 2 from Lemma 4. 

W.l.o.g. assume that u(<Si) > /'(iSa)- From upper bounds of circuit value we have 

A = r(5i)-t)(5 2 )^2 1+2 '. (16) 

Let R be a set of integers such that m € R iff m € [1, B\ and 

(v(Si) — v(S 2 )) mod m = 0. 

Note that if m ^ R then the algorithm outputs ‘v(<Si ^ u(<f> 2 )’. 

The set R cannot be large. At first we note that A is a multiple of lcm If r x. 
Suppose that #R > (l — e)B for e = (2f) _1 . We have e = (2f) _1 = (2 log 2 B)~ x < 
(2 In B )~ 1 . Thus, Lemma 4 implies 

A > lcm x£ i{ x > 2 n ( B / lnB ) = 

and we come to the contradiction with (16). This contradiction shows that #R < 
(1 — (2 £)~ 1 )B. Claim 2 follows from this bound. 


4 The proof of Theorem 2 

We prove Theorem 2 by adjusting an NC 1 -algorithm for integer comparison by Davida 
and Litow [5], So, we partially reproduce arguments from [4, 5], 

We start by some notation. Let pi = 3, • • • ,Pm be the first m odd primes. 

We denote the least nonnegative residue of x modulo n by [x] n : x = [x] n (mod n), 
0 ^ [x] n < n. For the rest of the section we set M = P 1 P 2 ■ ■ - Pm, Mi = M/pi , 
x i = [®]pi> Ci = [ x M i ] Pi . 

We denote the fractional part x — |_atj of x by {x}. 

By the Chinese remainder theorem any integer 0 ^ x < M is uniquely represented 
by Xi. For any integer x the following equality holds 

m 

Y J M i -tii = p(x)-M + [x\m. (17) 


Since Af;C < M, we have 0 ^ p(x) < m. The integer p(x) is called the rank of x 
with respect to {pi,... ,p m }. 

Dividing both sides of Eq. (17) by M we get 


E 


& 

Pi 


p{x) 


[x\m 

M ' 


(18) 


Eq. (18) gives a way to compute rank p(x) by approximating left-hand side. Let choose 
an integer h such that m/2 h < 1/4. For each i we take a (2 _ft )-approximationof £i/pi, 
i.e. 


fj_ < §i < s i + 1 

2 h ^ pi 2 h 


Sii h £ , 


where s* 



(19) 
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(20) 


Let a(x) = J2iL i s i2 h ■ Then, by summation of Eqs. (19) for all i we get 

{cr(x)} = p(x) — l&(x)\ + ~ a anc * 0 < a < 1/4. 

Now we reproduce Lemma 2.3 from [4], 

Proposition 1. The following holds. 

1. If{a(x)} ^ 3/4 then p(x) = [a(x)\. 

2. 7/1/4 < [*] m/M ^ 3/4 then |cr(x)} ^ 3/4 and p(x) = [a(x)\. 

3. If[x] M /M > 1/2 then p(x) = [a(x)\. 

All these facts are easily derived from Eq. (20). Using them we obtain an analogue 
of Lemma 2.5 from [4], 

Lemma 5. Let k*{x) = min(£; : k ^ 0 and {o([2 k x]M)} f 3/4). For any x the 
following holds. 

• k*{x) < oo. 

• k*(x) = 2P ol y( m ). 

• Suppose that k*(x) > 0. Then [x]m/M < 1/2 iff [[2 k x \m ]2 = 0. 

Proof. For x such that 1/4 ^ [ x\m/M f 3/4 we get k*(x) = 0 from Proposition 1.2. 
So, for this case the lemma is trivial. 

Suppose that [)x"j i\j /M *0 1/4 and k*(x) > 0. Choose k ^ 1 such that 


1 Mm 1 

2^+2 M 2 k+1 


( 21 ) 


Applying Proposition 1.2 we conclude that {cr([2 fc a;]M)} ^ 3/4. Therefore k*(x) ^ 
k < oo. The prime number theorem (15) has an equivalent form 


Pk ~ k In k. 


It gives a bound M = 2P°^ m \ Thus k*(x) = 2P ol y( m ). 

For k ^ k*(x) ^ k we use the upper bound of (21) to show that 2 k [x\M < M. It 
implies that [2 k ^[x]m]m = 2 fc ^[x ]m is even. 

Now suppose that [x]m/M > 3/4 and k*( x) > 0. Let y = M — [x]m = [v]m- 
So, [y\ M /M < 1/4 and [2 k y] M = [~2 k x] M = M - {2 k x\ M . Repeating the above 
argument we choose k ^ 1 such that 


1 < Mm < 1 

2fc+2 ^ M 2 k+1 


( 22 ) 


Since [2 k x] M = M-[2 k y\ M , we conclude that k* (x) = 2?°W"*) and [2^^[x] m ]m = 
M — 2 k * ( - x ' > [y] M is odd. □ 
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The algorithm described below compares integers represented by circuits over the 
division-free basis {+, —, •}. As was explained above, in this case comparing values of 
two circuits is reduced to comparing a circuit value with 0, i.e. to computing the sign 
of the circuit value. 

Let S be a circuit of size £. As it was mentioned above, |w(<S)| < 2 2 . There are 
to = 7t( 2 2<! ) — 1 primes from 3 to 2 2e . Due to the prime number theorem 

2 1M < to < 2 2£ . 

By the Chinese remainder theorem the residues modulo these primes represent uniquely 
all circuit values for circuits of size £. Moreover, |u(<S)| < M/A in these settings. 

The algorithm to check v(S) > 0. 

Input: S — a circuit of size l over the division-free basis. 

Step 1. Let h = 21. 

Step 2. Compute k*(v{S)). 

Step 3. If k*(v(S)) > 0 go to Step 6. 

Step 4. Compute ai := [u(«S)] 2) a 2 := [[u(<S)]m] 2 - 

Step 5. If ai = a 2 then output ‘yes’ else output ‘no’. 

Step 6. Compute b := [[2 fc *^( s ^t)(5)]Ar] 2 . 

Step 7. If b = 0 then output ‘yes’ else output ‘no’. 

Claim 1. The algorithm always gives a correct answer. 

Proof. Suppose that the algorithm finishes at Step 5. If w(6>) > Othentj(iS) = [u(«S)]m 
which implies [u(S)] 2 = [[u(<S)]m] 2 - If v(«S) < 0 then [u(<S)] 2 + [[u(«S)]m ]2 = [Mh = 
1. In both cases the algorithm gives a correct answer. 

Suppose that the algorithm finishes at Step 7. Note that in this case k*(v(S)) > 0. 
So, Lemma 5 can be applied. The algorithm gives a correct answer due to the following 
observation: x > 0 iff [x]m < M/2 assuming 0 ^ x < M/2. □ 

Claim 2. The algorithm runs in polynomial memory. 

Proof. Computing [u(«S)] 2 takes a polynomial time as mentioned above. The algorithm 
computes [[2 fe u(5)]M]2 by using Eq. (17) modulo 2 

m 

= p(x) + [x\ M (mod 2). (23) 

2=1 

To compute [x\m mod 2 we need to compute f, and p(x). Computing 2 k v{S) mod n 
is possible in polynomial time. To compute Mf 1 mod p, one can take all primes from 
3 to 2 21 one by one, compute an inverse residue modulo p, and multiply it by the current 
product. This as well as summation of can be done in polynomial memory. 
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For each i an (2 _2 ^)-approximation to <(, /p, is computed in polynomial time. Sum¬ 
mation of these approximations is directly implemented in polynomial memory. So, 
it is possible to compute [cr(2 fe z;(<S))J and {< 7 ( 2 * 11 ( 15 ))} in polynomial memory (for 
k = 2 poly W). Starting from k = 0 and incrementing k until k*(v(S)) is found may 
be implemented in polynomial memory. Having all these data, it is easy to compute 
p(2 fc *( w ( s ))u(<S)) and [[2*X 5 ))r;(,S)] M ] 2 . □ 

5 Open questions 

The main open question here is the complexity of SDFP. We suggest the problem AC 
as a ‘lowerbound’ for SDFP. If it is hard then SDFP is hard too. The maximal re¬ 
sult in this direction would be PSPACE-hardness of SDFP. Thus from the complexity 
viewpoint it is very interesting to put SDFP in PSPACEif possible. 

Complexity of AC^ itself is the next question. To be a good ‘lowerbound’ it should 
be hard. But up to our current knowledge nothing prohibits inventing an efficient al¬ 
gorithm for its solution. Such an algorithm would be interesting by other reasons. It 
would justify using circuit representation for rationals instead of binary system in com¬ 
plexity theoretic questions when time is estimated up to a polynomial slowdown. It 
might simplify complexity analysis of numeric algorithms very much. So, any definite 
result about ACj> would lead to interesting conclusions. 

From algorithmic point of view some intrinsic difficulty of numerical algorithms is 
related to nonconstructive nature of real numbers. Though some computational mod¬ 
els, e.g. BSS model of computation over arbitrary field initially proposed by Smale [2], 
directly operate with real numbers and thus completely ignore the question of number 
representation. Presently the relation of BSS model to the common model of algorith¬ 
mic complexity is not well understood. In particular, in BSS model SDFP falls into the 
complexity class NPr — an analogue of the class NP. Reffering to our way of number 
representation it is natural to ask about the inclusion SDFP £ NP AC >. This question 
is also open. 

In the opposite direction, it may be interesting to represent numbers as optimal 
solutions of SDP. For example, the well known n\ conjecture by Shub and Smale [15] in 
our setting asks about monotone circuit complexity of n\. It follows from the reduction 
in Theorem 2 that monotone circuit complexity is lowerbounded by the dimension of a 
SDP representing n! as an optimal solution. Is this dimension polylogarithmic on n? 
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